Numbers of donors and acceptors from transport measurements in graphene 
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A method is suggested to separately determine the surface density of positively and negatively 
charged impurities that limit the mobility in a graphene monolayer. The method is based on the 
exact result for the transport cross-section, according to which the massless carriers are scattered 
more strongly when they are attracted to a charged impurity than when they are repelled from it. 
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The discovery of graphene^ opened up a number of re- 
search directions* 2 The outstanding properties of this ma- 
terial originate from its massless electron-hole symmetric 
carrier dispersion e(p) = ivp, where the Fermi velocity 
v ~ 10 6 m/si 2 -i^^ Recent transport measurements^^ 
show that the graphene mobility is approximately inde- 
pendent of the carrier density. It is believed that the 
mobility is limited by the Coulomb impurities in the sub- 
strate, as well as by the quenched ripples.— In general, 
the carrier density and the mobility are sensitive to the 
immediate environment, e.g. to the presence of the ad- 
sorbed molecules, or to the donors and acceptors in the 
substrate. This has prompted proposals for applications 
of graphene monolayers as gas sensors Models^ ' 10 : 11 
based on the Born approximation for scattering off the 
Coulomb impurities, suggest a natural way to estimate 
the surface impurity density n,. 

In a typical situation, the Coulomb impurity density 
rii = n\ + rr[ consists of both positively and negatively 
charged species, here referred to as donors and acceptors. 
In this work I suggest a means to determine the densi- 
ties nf separately from a dc conductivity as a function of 
the gate voltage. This may find applications in the sens- 
ing of adsorbed gas molecules, that can be either donors 
or acceptors* 7 - as well as in the characterization of the 
substrate surface. 

The proposed effect is based on the recent 
observation^ 2 - that the exact transport cross-section 
is fairly sensitive to whether the charge carrier is 
attracted to an impurity or is repelled from it. For 
instance, for a conduction electron the transport cross- 
section is greater when scattering off a donor than off 
an acceptor. This attraction-repulsion asymmetry is an 
unexpected feature for scattering in the 1/r-potential, 
and is specific to the "relativistic" quasiparticle disper- 
sion. Indeed, such an asymmetry is absent in the exact 
nonrelativistic scattering off a 1/r potential in two^ 3 - 
and thresh dimensions, which is why the proposed 
effect would not work for conventional semiconductors. 
Remarkably, the asymmetry is also absent for scattering 
both off local (neutral) imperfections^ and off quenched 
corrugations^ Thus the proposed method selectively 
senses the imbalance between the charged impurities. 

Consider the scattering off the Coulomb potential 

Ze 2 ^ a 



where the dimensionless impurity strength 
a = Ze 2 Jhv, el = 2e 2 /(e + l) 



(2) 



can be both positive (attraction) and negative (repul- 
sion). Here Z is the impurity valence, and e is the 
dielectric constant of a substrate. The vacuum value 
Q! o|z=i.e=i ~ 2.2 for v w 1 x 10 6 m/s, while for the SiC-2 
substrate, ao|z=i,e=3.9 ~ 0.9. The effects of interactions 
between carriers in graphene, estimated via the scale- 
invariant RPA screening^ further diminish ao, 

a -> a = a /£RPA, £rpa = 1 + (tt/2) x el/Tiv . (3) 

This reduces the impurity strength by the factor 
£rpa|e=3.9 ~ 2.4 for a monovalent impurity on a SiC>2 
substrate, yielding \a\ « 0.37. 

The scattering states for the problem |T]) are char- 
acterized by the angular momentum j = m + 5, m = 
0, ±1, ±2, ... . The scattering phase shifts 5j arei 2 . 
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a e = a sign e . (4) 



Here 7 = y/ j 2 — a 2 , and |a| in this formula is limited by 
I (the subcritical impurit y 12 : 17 ). The property Sj = 5_j 
(mod 7r) ensures the absence of backscattering. Since 
the phase shifts depend only on signe but not on the 
absolute value |e| of the quasiparticle energy, the exact 
transport cross-section 



A tr = C{a e ) x A e 



9 - . 

C=-J2 sin 2 (S j+1 - S 3 ) (5) 
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is proportional to the energy-dependent carrier wave- 
length A c = 2nhv/\e\. The dimensionless function C(a e ), 
which is the transport cross-section in the units of the 
carrier wavelength, is plotted in Fig. [1] 

The transport cross-section (|S|) is strongly asymmetric 
with respect to sign a x signe, i.e. the sign of the po- 
tential as seen by the carrier: A donor (a > 0) scatters 
conduction electrons (e > 0) more effectively than it scat- 
ters holes (e < 0). The attraction- repulsion asymm etry 
is commonplace [e.g. Problem 6 to Sec. 132 in Ref . Il4| . 
Physically, one may expect the particle to spend more 
time around an attractive potential center and thereby 
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FIG. 1: (Color online) Even, C s (a e ) (red), and odd, C a (a e ) 
(blue), parts of the exact transport cross-section C(a e ) 
[Eq. ([5])] (black), their ratio C a /C s (dash-dot), together with 
the polynomial fits (If 5[) (dashed). Also shown is the Born 
approximation C Born — iLq? (green dashed). Inset: First ten 
scattering phase shifts 5j(a e ), j = Note that 5i/2 be- 

haves very differently from the rest. The differences Sj+i — Sj 
decay as Including only 61/2 and <5 3 /2 already gives a 

good approximation (red dashed) to the exact sum C (black) . 



be more significantly deflected (this intuition, strictly 
speaking, applies to the massive particles). However, for 
the practically important 2D and 3D Coulomb scattering 
in a parabolic band, the corresponding exact solutions 
are somewhat exceptional in a sense that they lack such 
an asymmetry. Remarkably, for the "relativistic" car- 
rier dispersion, characteristic of graphene, this generally 
expected asymmetry is recovered. 

The result (O applies to the half-filled 7r-electron band. 
In this case the RPA screening <|3j> is scale-invariant, pre- 
serving the functional form of the potential |T]). At fi- 
nite carrier density kp/ir, corresponding to the Fermi 
momentum kp, the RPA screening is not scale-invariant 
anymore i^iiS This means that the screened potential has 
the form fl} only at distances shorter than the screen- 
ing length l s ~ Z/akp, and is cut-off at larger distances. 
Scattering in such an effective potential is not tractable 
exactly, and the momentum relaxation rate is practically 
calculated only in the Born approximation ^ 9 ' 10 ! 11 The 
latter becomes asymptotically exact for small \a\ <C 1 
at large distances r ~ l s 3> 1 /kp, but misses non- 
perturbative effects at distances shorter than l s ~ 1/kp 
when the potential is large, a ~ 1. 

To investigate the role of the asymmetry, below I ne- 
glect the additional screening at finite kp, and use the 
scale-invariant screening ([3]) by the filled valence band 
only. Such an approach may be justified noting that, 
while the divergent Coulomb phased ln2fc^r (same for 



all the angular momentum channels) comes from the dis- 
tances kpr 3> 1, the relative scattering phases Sj, that 
contribute to the momentum relaxation via the trans- 
port cross-section, accumulate on the scale lj = j/kpJ£ 
This way the lowest phase shifts Sj are practically ac- 
quired within the screening length l s > lj , while the rest 
are strongly reduced by the screening at finite kp . Fortu- 
nately, the dominant contribution to the cross-section ([3]) 
comesi 9 - from the channels with j ~ 1 (Fig. [T] inset), that 
are least damaged by the additional screening at finite 
kp. In particular, the lowest shift 81/2 is behaving quali- 
tatively differently from the rest, causing the pronounced 
attraction-repulsion asymmetry. Such an anomalous be- 
havior is to be expected, as this is the channel in which 
the criticality, \a\ — j, is first reached. 

To calculate the conductivity, consider for simplicity 
the case when all the donors have the same valence and 
are characterized by the strength +a, while all the ac- 
ceptors have the opposite valence and the strength —a: 
hereon a > 0. Let /±(p) be the distribution functions for 
the particles and holes correspondingly. With the elec- 
tron interaction effects approximated via the screening 
([5]), the kinetic equations for each carrier type decouple, 

- df+ Sf+ * df- 6f- 

-et ■ ^— = ; — r , (6) 



dp T+(e p ) ' dp r_(cp) ' 

where e < is the electron charge, and 



v [n+A+ + n t A t 



[n+A tI . + n t A+] 



(7) 



Here A^ = A c x C(±a) are the right and left parts of the 
cross-section in Fig. [I] The rates ([7]) have straightforward 
meaning: The conduction electrons scatter off the donors 
with the enhanced cross-section Aj^ and off the acceptors 
with the reduced cross-section A^, while for the holes the 
situation is reversed. The Coulomb transport times ([7]) 
are proportional to the quasiparticle energy, 



-(e) = 
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2irhv 2 n+C(+a) + njC{-a) 
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(8) 



2irhv 2 n+C(-a) + n,-C{+a) ' 
From the kinetic equations (|6|) find the deviations 

5f ± = ±e(£p)vr±(e) x [-dj^] , p = p/p (9) 

of the distribution functions from the equilibrium Fermi 
distribution f±\e) — l/[e^ T>i ^ T + 1], where /i is chem- 
ical potential relative to the half-filled tt band. The re- 
sulting electric current [here Nf = 2 sp - ln x 2 va n oy = 4 
independent polarizations] 



eJ = eNfV 
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corresponds to the dc conductivity 
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Last, I express the conductivity (| 10[) via the carrier num- 
ber densities. Using integration by parts, the electron 
and hole densities n^ 1 can be cast in a similar form, 
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Combining Eqs. (JSJ), (flO|) . and (jTTjl . the conductivity 
(e 2 //i)n+ (e 2 //i)n 



+ nrC(-a) nttf(-a) + n;C(a) 



(11) 
(12) 



The asymmetry of the cross-section C(a) translates 
into the asymmetry in the dependence of the conductivity 
(|12p on the net carrier density n = n + + n~ . Introducing 
the charge carrier imbalance Sn — n + — n~ , as well as the 
symmetric and antisymmetric parts of the cross-section 

C s , a (a) = i[C(a)±C(-a)} ) (13) 

one can represent the conductivity (I12[) in the form 



a = 



n 1 
x — 



h niC s 



- (6n/n)c 
1-c 2 



/ \ Ca n f 

c(a) = J7- x — 



(14) 

Here rii = nf + n t is the total surface density of charged 
impurities. Thus measuring the conductivity asymme- 
try for, say, opposite gate voltages relative to the charge 
neutrality point, such that Sn ~ ±n, one can use fit- 
ting to Eq. (|14p in order to determine the concentra- 
tions nf of donors and acceptors separately, since the 
ratio c = (cr5 n= _ n -cr 5 „ =n )/(o'5 n= _ n +cr 5n=n ). Including 
the effect of the corrugations (ripples) 6 simply shifts the 
symmetric part niC s — ► niC s + C r in Eq. (|14[) . where 



a r = (e 2 /h)n/C r is the conductivity due to the rip- 
ples only. This still allows one to determine the donor- 
acceptor imbalance nf — nJ~ . To find nf + n~ one needs 
to independently gauge the ripple parameter C r by vary- 
ing the number of the charged impurities. As expected, 
the Born approximation C s —> C Born — ira 2 /2 has no 
asymmetry (C a — > 0), thus it entirely misses the effect. 

Fig- [H shows the odd and even components (fT3"|) of the 
exact cross-section, together with their polynomial fits 

C a w 1.68 a 3 + 7.69 a 5 -21.4 a 7 , (15) 
C s w 1.54 a 2 + 2.33 a 4 , C a /C s w 1.37 a - 1.06 a 3 . 

For |Z| = 1 impurities in SiC>2 (a ~ 0.37) the asym- 
metry C a /C s w 0.46 should be fairly pronounced. 20 Re- 
cent dc measurements of graphene monolayers on Si02 
substrate 3 show a noticeable asymmetry in the conduc- 
tivity a(n). The proposed explanation is that the ac- 
ceptors prevail over donors. This provides an alternative 
to the suggestion 10 that the asymmetry is caused by the 
gate-induced displacements of impurities. Further exper- 
iments with controlled numbers of adsorbed donors and 
acceptors may distinguish between the two scenarios. 

To conclude, a novel method for characterizing the 
electrostatic environment of a graphene sample is pro- 
posed. The method is based on the attraction-repulsion 
asymmetry of the exact transport cross-section for scat- 
tering off the charged impurities, and could potentially 
allow to separate the scattering effects of the donors and 
of the acceptors from those of the microscopic ripples. 
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